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CALCULATION OF SUPERSONIC FLOWS AT LARGE DISTANCES 
FROM SLENDER LIFTING BODIES* 

By Michael Schorling** 

Langley Research Center 

SUMMARY 

The exact gasdynamical equations are solved far from the axis of a slender body 
flying at an angle of incidence. The boundary conditions are obtained from slender -body 
theory. The solution allows one to make a prediction of the strength of the sonic boom 
as well as the position of the shock waves. 

INTRODUCTION 

Several theories have been developed for calculating sonic -boom ground pressure 
patterns and signatures. Most of these theories, such as those of Lansing (ref. 1), Hayes, 
Haefeli, and Kulsrud (ref. 2), and Friedman, Kane, and Sigalla (ref. 3), are based on geo- 
metric acoustics. These theories have great advantages in calculating sonic booms for 
unsteady flight conditions in an inhomogeneous atmosphere. Although corrections have 
been made by introducing nonlinear distortions, these theories are essentially linear in 
natvu:e and are not adequate for handling such nonlinear effects as shock focusing or cut- 
off Mach number. On the other hand, Whitham's theory (ref. 4) for the supersonic flow 
about noninclined axisymmetric bodies is based on a modified theory of characteristics. 
The results are obtained as solutions of the nonlinear differential equations of supersonic 
flow. 

The purpose of the present paper is to develop a theory to determine the supersonic 
flow in the far field of a lifting body with a nearly circular cross section. The nonlinear 
differential equations for supersonic flow are systematically expanded in a perturbation 
series. The resulting system of differential equations can be easily integrated for each 
order of magnitude. 


*The material presented herein is based on a thesis entitled "Berechnimg von Uber- 
schallstromungen weit entfernt von schlanken Flugkorpern mit Hilfe der Charakteristiken- 
theorie nach der P.L.K. Methode" submitted in partial fulfillment of the requirements for 
the degree of Doctor in Engineering, Rheinisch-Westfalischen Technischen Hochschule, 
i^chen, Germany, December 1969. 

**NRC-NASA Resident Research Associate. 


The flow of a perfect, nonviscous gas far away from a body flying at an angle of 
attack with a Mach number greater than 1 is considered. The body must be tapered, 
slender, and have a nearly circular cross section. Discontinuities in slope are allowed: 
however, abrupt changes in cross-sectional area (e.g., by a stabilizer) are not permitted. 
Small amplitude oscillations aroimd an axis normal to the direction of flight and through 
the center of gravity, as well as effects of an inhomogeneous atmosphere, can be included 
but are not considered in this paper. The extension to a stratified atmosphere is straight- 
forward. The extension woiQd allow an investigation of the nonlinear effects at the caus- 
tic. The flow is assumed to be steady, homentropic, and homoenergetic. There are no 
heat or mass sources in the flow. These conditions imply that the flow is irrotational. 


SYMBOLS 


ci = 1 - 


cross-sectional area of body, ttK^ 
speed of sound 

pressure coefficient 
n + 1 1 


„ _ o n + 1 sin fi 
C2 - -2 -n 


( 0 ) 




°r = |f 

r» , _ 9x 9 9x 9 _ Q . V . flu 
e^ = grad x^ 

e2 = grad x2 ) Cartesian base vectors 
e^ = grad x^ 


*m 


fimction of source distribution 


^m fimction of integration 
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H total enthalpy 

h specific enthalpy 

K(x, \l/) radius of body 

M Mach number 

m Fourier number 

m generatrix of Monge (or Mach) cone 

n ntunber of degrees of freedom of the gas 

n generatrix of wave normal cone 

n^o generatrix of wave normal cone in undisturbed flow 

n^ generatrix of wave normal cone in disturbed flow 

ng normal to shock surface 

p pressure 

R distance away from the body at which the boimdary condition is satisfied 

r cylindrical variable (radial distance) 

s entropy 

"s vector of the space 

T temperature 

t time 

t vector of the space E 

Uoo velocity of undisturbed flow field 
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w 


w 




X 


Xs 

xl 

x2 

x3 



magnitude of the speed, |wl 
vector of the speed 

components of the vector of the speed in a locally dependent base (i = 1,2,3) 
abscissa of characteristic surface 
abscissa of shock front 


= r cos xl/ > 
= r sin U/ 

J 


coordinates in the space 


= grad ^ 

= grad r = s’ ) vectors in the space 
= grad i// = 


z integration variable 

a ai^le of attack 

^ Prandtl factor, - 1 

y ratio of specific heats 

?7 Prandtl transformation, ^r 

Q angle of inclination of the streamlines 

X, Kj, /C2> factors of proportionality 

p. Mach angle 

I characteristic variable 

p density 
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II 


CT angle between shock front and undisturbed flow 

= iw,a,(p +\p) (i = 1,2,3) 

<p azimuth angle 

((> potential function 

\]/ cylindrical variable (azimuth angle) 

0 ) Prandtl-Meyer angle 

Subscripts: 


oo state of \mdisturbed flow 

* critical values 

Harmonic (or Fo\irier) numbers are denoted by numbers in the subscript position 
or by the subscript m. The order of magnitude is indicated by positive numbers in 
parentheses in the superscript position or by the superscript (j). The Greek letters p, 
V, a, and t used as subscripts and superscripts denote the Einstein summation con- 
vention. Single, double, and triple primes denote first, second, and third derivatives, 
respectively. A comma preceding a subscript denotes differentiation with respect to that 
subscript. 


ANALYSIS 

The exact nonlinear system of partial differential equations for supersonic flow 
is solved for large distances by using a perturbation method developed by Poincare, 
Lighthill, and Kuo (see ref. 5). 

The basic differential equations of the problem - the equations of continuity, irrota- 
tionality, and the condition of the characteristic surface — are expressed in invariant 
form. At first, this set of equations is expressed in cylindrical coordinates ix,r,\p). It 
is advantageous to express the velocity vector w in a particular locally dependent basis 
and to use the characteristic variable | in place of the independent variable x. Con- 
sequently, the differential equations are transformed to |, r, and \f/ as independent 
variables and to the magnitude w, the angle of inclination Q, and the azimuth angle cp 
of the velocity vector as dependent variables. The nonlinear equations obtained by this 
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procedure are split into systems of quasi-linear differential equations of different orders 
of magnitude by expanding each unknown function in a perturbation series. These sys- 
tems of equations can be integrated easily in each order of magnitude. The integrations 
lead to arbitrary fimctions of integration which are determined by the boimdary condition 
calculated from slender -body theory. All physical properties of the far flow field can be 
calculated by this method. By using the three-dimensional analogy of Pfriem’s formula 
(see ref. 6), the front shock can be calculated as well. Pfriem’s formula for two dimen- 
sions states that for weak shocks the shock wave bisects the angle between the character- 
istics before and behind the shock. 

Formulation of Basic Equations 

The following equations, written in their invariant forms, describe the problem. 
Equation (1) gives the condition of irrotationality; equation (2), the condition of continuity; 
and equation (6), the condition of the characteristic surface. 

Condition of irrotationality. - The condition of irrotationality is 


curl w » 0 (1) 

Condition of continuity . - The condition of continuity is given as 

^+pdivw = 0 (2) 

By the second law of thermodynamics ^T ds = dh - ^ dp^, the conditions of homentropic 
and homoenergetic flow are 


and 


H = h + i w^ 

The condition of continuity can be rewritten as 


a2 div w - ww • grad w = 0 


in which 


2 = ^ 
dp 


is the definition of the speed of sound. 
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Condition of characteristic surface. - Accordii^ to the theory of characteristics the 
generatrix of the wave normal cone n is normal to the generatrix of the Monge cone m. 
The quantity 4 denotes the variable along the generatrix n. Thus, the generatrix n 
is proportional to grad | ; that is, in oc grad ? or, with X used as the factor of 
proportionality, 

Xn = grad | (3) 

The speed of sound a is related to the velocity w by the relation 

w • (-"S) = a (4) 

Vs n is a imit vector (if . = 1), equation (3) yields 

X^ = grad ^ • grad | (5) 

ubstituting from equations (4) and (5) into equation (3) gives the equation of the charac- 
jristic surface in invariant form as 

a^grad 4 • grad ^ = (w • grad (6) 

The interrelation between the Monge cone and the wave normal cone is shown in fig- 
ire 1. The Mach angle is denoted by p. 

System of equations. - Within the system of basic equations composed of equa- 
ions (1), (2), and (6), there are five scalar differential equations for determining the 
mknowns | and w. Since the flow is homentropic, there exists the relation p cc py 
oetween pressure and density where y is the ratio of specific heats. Hence, 

h = na^ + Constant 

where n = denotes the number of degrees of freedom of the gas. Combining this 

equation with the equation for homoenergetic flow gives 

na^ + w^ = Constant 

Thus, the speed of sound a is a function of the velocity w. 

If one thinks of equation (6) as the defining equation for then there are four other 
equations (eqs. (1) and (2)) for determining the vector of speed w. Thoi^h the three 
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equations curl w = 0 are dependent on one another, they must be used together with the 
condition of continuity (eq. (2)). Alone, curl w = 0 gives the solution w = grad as 
curl grad (f> = 0. 


Change to a Characteristic Coordinate System 

After having formulated the system of equations in an invariant form, one must 
select a coordinate system in which to carry out the solution. A plotted point P has 
the coordinates 


x2 = r cos ij/ 


and 


X 


3 


= r sin 




as indicated in figure 2. 

The imit vectors in the space are 

"e^ = grad 

e^ = grad x^ ) 

e^ = grad x^ 


(7) 


as shown in figure 2. From equation (7) one 
grad r = s^ = cos + sin 
grad j// = iT= ^(-sin + cos i 

As shown in fig\are 3, e^,^,T is also a system of orthogonal vectors. It is useful 
to express the velocity vector w in a locally dependent basis and to define a new set of 
velocity components in this basis. The space E^(e^,s,T) just introduced will be used. 
According to figure 3 the velocity vector can be related to its magnitude w, its ar^le of 
inclination 3, and its azimuth angle cp + x}/. So the velocity vector is written as 


calculates 
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(9) 


w = wjcos 


+ sin a(cos <ps + sin 



In the space E^(el,'s»^) the new components of velocity Wj, Wg, and Wg are defined 
as 



(i= 1,2,3) (10) 


where 


oi = (w,3,<p + if/) 


(i = 1,2,3) 


So far, the unknown fvmctions have been expressed in terms of cylindrical coordi- 
nates; that is, 

w = w(x,r,^) 


and 


? = l(x,r,j^) 

Since disttmbances sent out by a body travel along characteristic surfaces, it is preferable 
to introduce a "characteristic variable" $ along the generatrix "H. Hence, the indepen- 
dent variable x is replaced by the formerly dependent variable This step follows 
the Poincare -LighthUl -Kuo (P.L.K.) method (see ref. 5). 

The new independent variable defines those characteristic surfaces propagated from 
the body. Disturbances traveling along the other family of characteristic surfaces are of 
lower intensity. For this reason, and also because the coordinate r seems to be phys- 
ically much more reasonable, it is convenient to retain the variable r and avoid intro- 
ducing a second characteristic variable. Now 

w = w(|,r,j//) 


and 


x = x(4,r,)^) 
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Final Form of Equations 

The base vectors y^, and y^ for the coordinate system |,r,i// are defined 
by the equations 


yl = grad | 

y2 = grad r = S’ S 

y3 = grad \f/ = 

Although "y^ is a unit vector, y^ and y^ are not. 

In the space Y^(yl,’y2,y3), 

el = grad xl = e^ = grad = x^yl + Xj.’y^ + x^y3 


( 11 ) 


( 12 ) 


Here and in the following equations, the Einstein summation convention is used. The 
Greek subscripts and superscripts /i, ct, and t indicate this summation convention. 

From equation (12) one determines that 




■) 


Analogously, the operator grad is determined in the space as 


(13) 


X| grad = elD^ +"sDr +^TD;^ 


(14) 


where D^, Dp, and are abbreviations for the operators defined by the equations 

D| = Dp = x^ 8p - Xp 8| = d^jy - xp 8| 


Furthermore, 


curl w = e^ X -^1^ = grad x Wq- 
8u^ 8x^ ^ 




div w = e”^ • = grad • w^. 

8u^ 8x' ^ 


(15) 


10 



By usli^ equations (15), equations (1), (2), and (6) can be written in terms of the velocity 
components, wj, W2, and W3 and the independent variables r, and \l/. The 
results are 


w j X grad w + W2 x grad ^ + W3 x grad(<^ +il/) = 0 

-cot^H wj • grad w + W2 • grad a + W3 • grad(<p +)//) = 0 > (16) 

(w • grad 1)^ = a^grad | • grad ^ ^ 

To simplify eqxiations (16) the Prandtl-Meyer fvuiction do) = cot ju ~ is intro- 
duced. Here, o) denotes the Prandtl-Meyer angle and p denotes the Mach angle. By 
using equation (10) the vectors wj can be expressed in terms of the orthogonal vectors 
e^, "s, and T. The five scalar differential equations obtained by this procedure are 

Condition of irrotationality: 

In direction of e^, 

tan p sin 0 sin i tan p sin 0 cos cpDipo) + cos 0 sin <pDrO 

- ^ cos 0 cos (pD^pd + sin 0 cos cpDr(p + ^ sin 0 sin <pD^<p + A sin 3 sin cpx^ = 0 

(17a) 


In direction of s, 

tan p sin S sin <pD|W - — tan p cos + cos 0 sin <pT>^Q 

+ A sin + sin 0 cos cpD^cp = 0 (17b) 

In direction of T, 

tan p sin d cos cpD^u} - tan p cos ^Dr(o + cos 5 cos 

+ sin SDi-a - sin S sin cpD^cp = 0 (17c) 
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Condition of continuity: 


-cot n 


cos 5D|W + sin 5^cos cpBj.u) + p sin - sin + cos a^cos cpDj.^ 


+ i sin (pl>^ - sin O^sin (pl>j.cp - ^ cos cpD^p^J + i sin 3 cos <px^ = 0 (17d) 


Condition of the characteristic surface: 

2 

^cos ^ - sin 3 cos ^Xj. - sin 3 sin ^ i x^) = sin^/i ^1 + Xj.^ + x^^j (17e) 

These five equations are the starting point for a perturbation calculation. 

Calculation of Perturbation 

The system of nonlinear partial differential equations (eqs. (17)) is solved by using 
a perturbation method. The following assumptions are made for the unknown fxmctions: 


CO = + . . . 

X = x^®) + X^^) + X^^) + . . . 

cp = (^( 0 ) + (^( 1 ) cp{^) + . . . 




(18) 


The superscript in parentheses denotes the order of magnitude. It is assumed that for 
large distances r the zero -order term is greater than the first-order term, the first- 
order term is greater than the second-order term, and so forth. These series do not 
contain any assumption about the value of the orders of magnitude. The terms in the per- 
turbation series can be calculated without making any such assumptions. 

Some prior knowledge about the far flow field to the zero order can be used to 
shorten the calculations. For large distances r. 
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o)(0) = Constant 


= 0 


= Constant 


= Constant 


0T* rr 




dj. « dt 


9j» 3|i 


but Xp x^ ~ X|. The assumption = 0 implies that the abscissa e^ = grad x^ is 
in the direction of the undisturbed flow. 

It is advantageous to develop each of the terms in the perturbation series in a 
Fourier series in ij/. Convenient forms for these developments are 


oo 


(i S 1) 


(20a) 


where B stands for o), aJid so forth, and 




(j ^1) 


(20b) 


Naturally, the system of equations developed in Fourier series must also be satis- 
fied for each Fourier number m. By developing each of the five scalar differential equa- 
tions (eqs. (17)) into Fourier series and suitably combining equations, one gets equations 
for the different orders up to the second order: 


"mj? tsSJl - 0 


,<1) . 0 


.Jl) .o(l) +i3(l)_o 

^m,r ^ 9m,r + r ^m “ ^ 


(21a) 


(21b) 


(21c) 





tail 

+^m,Vtan 

+^S^) = 

^m,r + ^mlr + F = F ~ f <i^m-i/^i^^ 

sin2/x^0)xj^)^ = + sin^pt^^^cos /^^°^x{^!-i^,r^ijr 


(21e) 

(21f) 

(21g) 

(21h) 

(21i) 


In these equations, the Fourier number m takes integral values in the range -<» s m ^ ®o^ 
and the summation convention is indicated by the Greek letter v, with -oo ^ 2 oo. in 

fact, for axisymmetrical body shapes one can restrict v to 0 ^ v ^2. Within this 
restriction the desired functions are sufficiently approximated. 


Zero-, First-, and Second-Order Perturbation Solutions 

In deriving the system of equations (21), the results of zero -order calculation have 
already been used; that is, 

x^®^ = cot 


and 


cp(fi) = 0 

Because of the special structure of equations (21), one cannot obtain any information about 
In linearized two-dimensional flow the characteristics are straight lines parallel to 
each other and obey the law 


X = I + r cot pi 

From this equation one estimates that in three-dimensional flow 

x(0) = 1 


(22) 
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By using the relationship between the Mach angle fx and the Prandtl -Meyer 
angle w in homoenergetic flow, one gets 


w = r cot M ^ f cot p, = M + (n + l)tan-l[-J^^=^'\ - | (23) 

-^w ^ ‘^M* \^n + Ij 2 

From a Taylor series for p in powers of o), one obtains the following expression' 
for IX, which is correct to second order: 

jx = ix^^^ + i 02 ^^^^ (24) 

Here Cj and C 2 are constants defined by 


N = 1 - " •*• ^ 
'1 n 


COS^jU.^®) 


and 


= 2 n + 1 sin M 


( 0 ) 


n 


cos^ix^^^ 


With the use of equation (24), this system (eqs. (21)) can be integrated in a system- 
atic step-by-step procedime. The details of the integration are given in reference 7. 

The first- and second-order solutions are as follows: 

^m = 

^m = -GJi^(4)r-l/2 

IxS^ = -ciG^\|)r-l/2 ^ (25a 

x^^sin^ju,^®^ = -2^1 - Cj^G^\|)r^/2 




= -cot 

+ (| - ^)tan jii ^° V - 3/2 + Gj^\^)r-l/2 

= cot M^°^Gjj!.p(^)G|f^(4)r-l 


^ 


+ li+Si^ 
\8 2 


)tan M(0>r-3/2 jj G{l)U)d^ - Gg)(|)r-l/2 

^i{2) = cj cot +1 C2 |g^^(^)] r'^ 

+ cj(^ + ^)tan ;x^°^r-3/2 Gi^\|)d| - ciGj^\|)r-l/2 

^{2)sin2|Li^0) = ^1 - ci)(2 + cjcot 

- f C2 [gS’( 5^ V r - 2(1 - Ci)og>({)rl/2 


+ tan f -(l - ci)m2 + 1 (3 + cj) 


,- 1/2 




(25b) 


The terms G^^(^) are functions of integration with respect to r which are to be cal- 
culated from boundary conditions. 


Coordinate of Shock Waves 

In order to calculate the coordinates of the shock wave in the space Y^, Pfriem's 
formula is extended to this problem. Pfriem’s formula of plane steady flow states that 
for weak shocks the shock wave bisects the angle between the characteristics before and 
behind the shock (ref. 7). One can extend the formula to the space y 3 by considering 
not the shock wave but the normals upon the shock surfaces. 
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The followii^ definitions are used: 

n^o normal to the generatrix of the Moi^e cone in the imdisturbed flow 

n^ normal to the generatrix of the Monge cone in the disturbed flow 

lig normal to the shock surface 

These vectors are shown in figure 4. 

The extension of Pfriem’s formula to space is given by 

Hoo • ns = “I • ^s 

Analogously to equation (3), set 

Hoo = Kj grad 

1i| = ^2 Srad | > 


(26) 


(27) 


\ = /C3 grad J 

with Kj, K2> 3-nd as factors of proportionality chosen so that n< 5 o, n^, and iig are 
unit vectors. 

Combining equations (26) and (27) yields 


sin^ 




(28) 


Equation (28) must be developed in the same way as equations (17) to obtain a set analo- 
gous to equations (21). Up to the second order, this set is 

xs,r<°^ = cot p(0) 


(1) - lx (1) - - 1 (3(1) + n(D) 

3,r - 2 *r " . 2 (0)'^ ^ ^ 

2 sin^ja' ' 

g^j.( 2 )sin^p(®) = - 4(9^^^ + + sin^jUL (O)cos p(®)^Xg 


(29) 
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These equations can be integrated to give 
= r cot ^ 

Xsm^sin = -(l - Ci)GjJ\^)rl/2 

Xgfsin2M<‘'>=-(l-ci)G<%rV2 

+ + Ci)cot ^ 1“ r 

For general three-dimensional flow, the solutions are valid up to first order; for axisym- 
metrical flow, up to second order. 

Satisfying the Boundary Condition 

In equations (25) and (30) the functions of integration G^\|) were introduced. The 
functions G^^d) and G^\|) must be calculated from the boundary conditions. Unfor- 
tunately, the boundary conditions are formulated at the body itself thoi^h system (21) has 
been integrated for large distances r away from the body. In this paper, the boundary 
condition will be obtained from slender -body theory, a linearization of the singularity 
method. Use of this linearization is another approximation. But slender -body theory has 
the advantage of short, compact solutions. As system (21) has been integrated for large 
distances r, the solutions cannot be continued to the body itself to satisfy the boimdary 
condition there. One must satisfy the boimdary condition at some distance R away from 
the body. For this purpose, the following procedure is used. 

It is assumed that the disturbances sent out by the body run along straight charac- 
teristics (zero -order theory). At the distance R the straight characteristics are 
matched with those characteristics which result from the exact theory and which are 
valid for large distances. A discussion of the matching distance R is presented sub- 
sequently in connection with an example calculation. 

From a physical standpoint, a proper way to calculate the functions of integration 
is by considering the angle of inclination of the streamline. The inclination may generally 
be expressed by 



tan 3 = 
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Here v denotes a vertical velocity component and Uoo + u, a horizontal velocity com- 
ponent. For small disturbances this equation can be approximated as 


a = 



^ 8 r 




dr) 


/ 3 <#) 


V 


(31) 


where /3 = ifM^ - 1 and 77 = /3r. The boundary condition is satisfied by usii^ the 
equation. 

Developed in perturbation and Fourier series, equation (31) can be rewritten as 




2 cos \l/ + 2 cos 2i// 

+ + <#> 77 f ^)2 cos yp + 


+ cos 2rf/ 




(32) 


By comparing equation (32) with equations (25), the function of integration can be 
expressed in terms of tp^j, which itself can be determined by slender -body theory. The 
comparison gives 




= ^^<Pr)Q^ + cot - |tan G^>^^4)d4 


g(2)(|) = P/ 3 (#>^( 2 ) + Jtan G^^^>(?)d4 


(33) 


2 

G^2)(|) = ^^<^^(2) ^ M^°^^_^^(?)Gj,^>(|)R-l/2 +^tan m<0)r-1 G<l)(|)d| 



In order to determine the quantity a modified equation given by von K^man and 
Moore (ref. 8 ) is considered: 


0 


oo 

L Jz=o 


^m 


(z) 

' ' m 


m=-oo 


Here fm(z) denotes source distribution functions (which are given subsequently) and 


sinh 


X = l/cosh^X - 1 = 


The derivative with respect to rj = /3r gives 



(34) 


Developing the integrals in equation (34) for slender bodies, or small rj, gives the fimc- 
tions fjn* Since only symmetrical bodies are considered, one can restrict m to 
0 = m ^ 2 and obtain 


= T'” *0 dz - f'” f, -£2 sA ^ 3-1^) 

f cosh 2X cosh X dz(e2i«/^ + e-2i»^) 
jQ ^ V sinh X ' 

The fimctions fjn(z) which follow from slender -body theory are 

^ fo(z)dz = - v)^oi^ - V) 


fl^"2)(x - 77 ) = - ^ Ki’(x - 77 ) - ^ 


(35) 


(36) 


- 77 ) = - K 2 ’(x - 77 ) 
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The angle of attack of the body a is measured counterclockwise, as shown in figure 5. 
The radius of the body is denoted as K(x,i//) and can be developed in a Fourier series 

oo 

as K(x,i//) = ^ Kjji(x)e^“’^. Since one supposes equations (36) to be introduced in 

m=-«> 

equation (35) and, after that, equation (35) developed in first and second order, one is able 
to calculate the functions of integration (eqs. (33)). This procedure shall now be demon- 
strated for a body of revolution. 


The Body of Revolution Flying at an Angle of Attack a 
For a body of revolution the radius is simply 


K(x,.//) = Kq(x) 

Therefore, the source distribution functions (eqs. (36)) are 

fo^"^^(x - r?) = -Ko(x - t?)Ko’(x - t?) ^ 

or \ 

fo(x - ^) = - ^ - V) 

2 

where A-ttK denotes the cross-sectional area of the body and 

or ) 

fl(x - ’l) = - ^ a^A”(x - 77 ) 


(37) 


(38) 


The source distribution functions with higher Fourier numbers are all equal to zero. 

Thus, for a lifting body of revolution, only the Foiarier numbers m = 0 and m = 1 have 
to be taken into account. By putting these nvimbers in the set of solutions given by equa- 
tions (25) and referring to equation (20a), it can be seen that 
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or because of the symmetry at i// = 0, 

^cos i}/ 

The solutions can then be determined. 

The solutions for S and x are written as 


cos i// = + 2 cos 


sin^jLi' ' ^ 


(39a) 


and 

^(2) _ ^ 25^^^cos ij / 




-cot /x'^^'r 


(O)r-l 


,(1) 


+ + 4 cos 


+ + 2 cos 


+ i tan jx(®)r 
o 


3 jj Gf)^\|)d| - 2 cos G5_^\|)d| 




= <||l - Cjj^2 + Cj]cot 


+ 4 cos 


1^2 


+ 2 cos i//(gS^)| 


) (39b) 


>ln r 


-2(l-Ci)rl/2(G(2). 2 cos xp'G^^ 


+ tanM<0)^^^r-l/2 


- tan - Cjjr"^/2 G^^\^)d| 
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Calc\ilation of the Functions of Integration 

For a lifting body of revolution, the functions of integration (eqs. (33)) can be writ- 
ten as 





+ cot 




I tan G^(|)d| 


and 

+ 2 cot +itan G^^^^d? 

1 vj 0 1 o 

Here, the quantities cj) are to be determined. This determination can be made, as 

”m 

mentioned previously, by using the slender -body theory (eq. (34)). By using the derived 
source distribution fimctions for a lifting body of revolution (eqs. (37) and (38)), equa- 
tion (34) can be written as 


_ 1 (X - z)A"(z) ,, 

^0 27rr?Jo ^(x - z )2 - r ?2 

> 

^ _ Q!)3 (x - z)^A"(z)dz 

’’I ^-^0 _rP‘ ^ 

or, because ^ + ?7 = x, as 

^ _ _1_ ('^ U +V - z)A"(z)dz 

27777 Jq + 2 t 7 - z 

. _ dig (g + 77 - z)^A"(z)dz 

27777^ *^0 l/g - z /g + 277 - z 


(40) 


(41) 


In calculating the functions of integration 77 means here 77 = ^R. 
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By integration by parts the singularity at z = | can be removed. The result is 


2H7} 


.21^I|±£^A"(z) 

^ + 2?7 - z 


1^1 + 2?7 - z 


Jo ^ i/r I ^ 


l/(? + 2 t; - z)^ 

After the development of the integrand in power series, the integration is possible. As 

a result of these operations, the first term of (j>„ vanishes. Since the values of rj 

'0 

considered herein are large (77 > 10) and the values of | are small (| = 0 to 1), it can 
be seen that the magnitude of the second term of is about 77 times larger than the 

third term of Therefore, the second term is of the first order and the third term 

is of the second order. One gets for pointed bodies (A(0) = A'(0) = 0): 


<t>. 




1 

\UU+v)l 


A1 


. 5A3 

271-77 

\lg+2n\ 


2^2 

sr 

4 ■ 

16|^ 

1 

'^(1 + 377 ) 

/a'. 

A1 

A2 

A3 . 

217-77 

\lig + 2 t 7)3 

r 

'2g 

8^2 

16^3 ^ • 


(42) 


Here and in equations (43) and (44), 

A’»M=27tK^ 
d^ ag 


A1 = A’l - A 

A2 * A’?2 _ 2Ag + 2 A(z)dz 


and 


A3 = A'|3 - 3A|2 + 6 A(z)z dz 

Jn 
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Analogously, the term is calculated to be 


Hence, 


, _ gg (g + 77 - z)^A"(z)dz 

217772 Jq - z + 277 - z 


. m_ \/F(l +v)^ Ia' . A1 ■ 3A2 . 5A3 . ) 

277772 + 2t7 2^2 8^3 16^4 ‘ * y 

* 2^(i «-2,)3 \ 8{2 la«3 ;j 

The functions of integration are determined as 


_ ifFd + T?) /a’ 

3 . A1 ^ 

2i7|/R (/I + 277^1 

■I- — 9 + 

2^2 

_ a^(?+77)^ /a 

’ + Al. 

2i7R^/^^ + 277 V 

2|2 


(43) 


(44a) 


and 


2i 7^^(| + 277)3 \ 8|2 16^3 y 


+ cot fjL 


(0) 






„^2) . «03e‘ *mi * 7 ^)/ ., . Al . A2 A3 . ^ 
4;7rS/2|/(S + 2,)3 \ 8«^ 16«’ ) 


+ 2 cot 




(44b) 
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The Sonic -Boom Signattire 

Finally the sonic -boom signature is calcxilated. The pressure coefficient is 


P-Poo 

-T - Cp 2<pj^ 

JL 'll ^ 

2 


^Poo'loo 


where is the disturbed velocity in dimensionless form in the direction of e^ and 
where the subscript <» characterizes the values of the free -stream flow. The disturbed 
velocity <^x 

<#>x = £ cos a - 1 

Hence, the pressure coefficient becomes 


Using the energy equation of homoenergetic flow 
naoo^ + Uoo^ = (n + l)a*^ 
and the definition of speed of sovmd 


2 _ n + 2 Po 


3-00 ~ 


one gets 


P P°o _ _ n + 2 


n 


cos 


i/n + Mo 

|~irrT 


2 

■ Mo 


( 45 ) 


dM, 


and is a function 


The term M* is given by the Prandtl-Meyer angle o) = \ cot ju. 

Jn M* 

of jx and o). The integration is possible because p can be expressed as a function 
of Mile. The integration gives 


M* = 


/co - jui + ^ \ 

„ (n + l)tan2 ^1 + 1 

2 \ n + 1 




tan2 


w - M + j 

I n + 1 


+ 1 
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where 


H = = sin"^ ^ C20>^^^ 

and where 

(t) = 0)^®^ + 0)^^^ + 

The term 0,(0) is given by 

0,(0) = n(0) + iSm tan-1 ^^^ ^ ^ - | 

and w(^^ and u)(^^ are given by equations (25). 

In equation (45), 

a = ^(1) + ^(2) 

and 

cos ^ = cos (^(^^ + ^(^^) 

The terms ^(^^ and ^(^^ are also given by equations (25). 

Comparison With Result of Reference 4 

For a body of revolution at zero angle of attack, the equation of the abscissa x is 


= r cot p(®^ + 1-2 — 


sin2p(0)cos2p(0) 


g(1)(|)|(F 


- 2 


n + 1 


sin^/x(0) cos^ju,^^^ 


G^q\^)\1¥ 


4n cos^/i(0^ - n - 1 ,^(1) 

4n cos^p(0)sin 


2n + 1 n + 1 


n 


^<’>(4^1 


^ cos3ju(0)sin3p(0) 
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By using n = — and - 1, the previous equation becomes 

V - 1 


X = r^Mo 


1 + « - (r - 1) - (y * « "7 

M J - 1 “ uj - 1 “ 


Mqo p(^)/ 


(m«2 _ ^ I 

^ > L G<«(«)dJr-l/2 


(m«, 2 . l)' 


(y + l)(y + 3) Moo 


6 


{mJ - 1)' 




In r 


In comparison, Whitham’s result (ref. 4, eq. (99)) is 


= t\ImJ - 1 + ^ - (y + 


1) ,, f(4)VF 

7Moo^ - 1 


== ^^-|-^[Moo^(y + 4) - (2y + 5^ |F(^)]^log r 

l) 

Though it is impossible to equate the terms of both solutions as different methods 
of calculations have been used, similarities in coefficients can be seen. Above all, the 
comparison shows that, in general, both solutions obey the same power series. But one 
must recognize that in the present paper no approximations have been made in obtaining 
the solutions of the nonlinear differential equations; however, approximations are made 
in evaluating the functions of integration. The slender -body theory has been used in both 
papers to calculate the functions of integration. An extension of Whitham's theory for 
slender lifting bodies is given in reference 9. 
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EXAMPLE CALCULATIONS 


The sonic -boom pressure signature has been calculated for different Mach numbers 
for two bodies of revolution. The cross-sectional areas of the two bodies are given in 
dimensionless form by 
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Aj = 0.04(x^ - 2x^ + x^)ir (0 = x = 1) 

and 

Ag = 0.0i(4x 2 - 4x^ + x4)it (0 £ x i i) 

The numerical evaluation is made by using equation (45). 

In figure 6 the fimctions of integration and are drawn for the two 

body shapes considered. As was mentioned previously and is illustrated in this figvire, 
the second-order integration function is much smaller than the first-order integration 
function; this fact provides a justification for restricting oneself to a first-order calcula- 
tion. Thus, the influence of the second-order calculation has been neglected in the exam- 
ple calculations because of the appropriate choice of the first- and second-order integra- 
tion functions (eqs. (44)). 

The boundary condition was satisfied at some distances R away from the body. 

As shown in figure 7, the function of integration G^^\|) is independent of the value of 
the matching parameter R over the range 3 i R = 50. For smaller values of R, the 
solutions broke down as the set of differential equations was solved for large distances r, 
whereas for larger values of R, the boundary condition was not accurately represented 
by slender -body theory. 

As shown in figure 8, the point of intersection of the Mach line of the undisturbed 
flow through the nose of the body and a parallel line to the abscissa x at the radial dis- 
tance r gives the origin of the plots of the pressure signatures shown in figures 9 to 11. 
In figure 9 different pressure signattares for one body of revolution are drawn. The angle 
of attack at, the azimuth angle \p, and the radial distance r where the pressure signa- 
ture is calculated are fixed. Here, the pressure signature is a function of the distance R 
at which the boundary condition is satisfied. It can be seen that the different values of the 
matching parameter R do not influence the shape of the pressxire signature. In all the 
sample calculations only the bow shock was calculated according to equations (30). 

Figures 10 and 11 show the influence of the factor a cos p and the radial dis- 
tance r on the pressure signature. Naturally, the pressure rise increases with 
decreasing distance r and increasing Mach nvunber. Fvtrthermore, the pressure rise 
for a lifting body is greater at the bottom side ()// = 180°) than at the top side (p = 0°). 

In figure 12 a comparison of the results computed by the present theory and those 
computed by Whitham's theory (ref. 4) is shown. The agreement between these theories 
is rather good. However, it must be emphasized that the present theory has advantages 
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such as simplicity of equations, handling of lifting bodies, and easy extension to more 
general problems as the change from a homogeneous to an inhomogeneous atmosphere. 

CONCLUDING REMARKS 

The exact nonlinear system of partial differential equations for supersonic flow was 
solved for large distances by using a perturbation method developed by Poincare, Lighthill, 
and Kuo. The unknown functions were expanded in a perturbation series. The system of 
differential equations were derived for each order of magnitude and the unknown functions 
obtained by step-by-step integration. The integration introduced arbitrary functions of 
integration which could be used to satisfy boundary conditions at the body. The boundary 
conditions were obtained as in previous studies of nonlifting bodies by using slender -body 
theory. However, axisymmetric lifting bodies with nearly circular cross sections were 
considered herein. In order to calculate the flow away from the body, the boundary con- 
ditions were satisfied not at the body itself, but at some distance R away from the 
body. The example calculations give reasonable sonic -boom signatures for arbitrarily 
chosen bodies of revolution. The results are in good agreement with the previous com- 
putations for nonlifting bodies. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 20, 1971. 
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Figure 4,- Description of the weak shock (Pfriem’s formula). 
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Figure 5.- The shock front and the characteristics. 
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(a) Aj = 0.04 (x 2 - 2x^ M= 1.1. (b) Aj = 0.04(x2 - 2x^ +x^)ir; M= 2.5. 

Figure 6. - Comparison of the functions of integration and R=20; acosi//=0.1, 









Dimensionless length x 


44. D 


Dimensionless length x 


(a) a cos ^ = 0; r = 50. 
Figure 10. - Influence of the distance 


(b) a cos i// = 0; r = 100. 

r on the pressure signature for A 2 = 0.01^4x2 - 4x^ + 
M=l.l; R=20. 
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